arXiv:1508.00770v2 [cond-matsoft] 14 Oct 2015 


Tailoring superelasticity of soft magnetic materials 

Peet Cremer.^'Eil Hartmut Lowen/ and Andreas M. MenzeP F^ 

Institut filr Theoretische Physik II: Weiche Materie, Heinrich-Heine-Universitdt Diisseldorf, D-40225 Diisseldorf, 
Germany 

(Dated: 15 October 2015) 

Embedding magnetic colloidal particles in an elastic polymer matrix leads to smart soft materials that can 
reversibly be addressed from outside by external magnetic fields. We discover a pronounced nonlinear su¬ 
perelastic stress-strain behavior of such materials using numerical simulations. This behavior results from a 
combination of two stress-induced mechanisms: a detachment mechanism of embedded particle aggregates as 
well as a reorientation mechanism of magnetic moments. The superelastic regime can be reversibly tuned or 
even be switched on and off by external magnetic fields and thus be tailored during operation. Similarities to 
the superelastic behavior of shape-memory alloys suggest analogous applications, with the additional benefit 
of reversible switchability and a higher biocompatibility of soft materials. 


The term “superelasticity” expresses the capability of 
certain materials to perform huge elastic deformations 
that are completely reversibleii^. It was initially intro¬ 
duced in the context of shape-memory alloys^Ti^. These 
metallic materials can perform large recoverable deforma¬ 
tions due to stress-induced phase transitions. A transi¬ 
tion to a more elongated lattice structure accommodates 
an externally imposed extension. Typically, this transi¬ 
tion shows up as a pronounced “plateau-like” regime on 
the corresponding stress-strain curve. On this plateau, 
the samples are heterogeneous with domains of already 
transitioned material. Then, only relatively small ad¬ 
ditional stress induces a huge additional deformation. 
Smart material properties are observed^”—: upon stress 
release, shape-memory alloys can reversibly find back to 
their initial state. They self-reliantly adapt their appear¬ 
ance to changed environmental conditions. 

In the present letter, we demonstrate that an anal¬ 
ogous phenomenological behavior can be realized for a 
very different class of materials, exploiting different un¬ 
derlying mechanisms. Moreover, we show that dur¬ 
ing operation the behavior can be reversibly tailored 
from outside by external magnetic fields. All of this is 
achieved by employing soft magnetic gels as working ma¬ 
terials: colloidal magnetic particles embedded in a possi¬ 
bly swollen elastic polymer matrix—. Similarly to mag¬ 
netic fluidsii“— , magnetic gels allow to reversibly ad¬ 
just their material properties by external magnetic fields. 
In this way, switching the elastic propertiesi^Ti^^ offers 
a route to construct readily tunable dampers^i or vi¬ 
bration absorbers^^, while the possibility to switch the 
shap e allows application as soft actuators^Sri^. 

Here, we show that magnetic gels due to the inter¬ 
play between magnetic and elastic interactions likewise 
feature superelastic behavior: it is enabled by a de¬ 
tachment mechanism of embedded magnetic particle ag¬ 
gregates and by a reorientation mechanism of magnetic 
moments. Both mechanisms are stress-induced and re- 
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spond to external magnetic fields. Therefore, superelas¬ 
ticity can be switched on and off, and also its magni¬ 
tude and position on the stress-strain curve can be re¬ 
versibly tailored during operation, as has been realized 
for some special metallic components^i^. The supere¬ 
lastic behavior in our case covers a significantly larger 
strain regime. Furthermore, soft gel-like materials gener¬ 
ally provide a larger deformability and higher degree of 
biocompatibility^-— than metallic alloys. This becomes 
particularly important for medical applications^. There 
has been significant effort to transfer the properties of 
shape-memory alloys to soft material o^^i^^i^^dQ . Here we 
report on reversibly tailoring superelastic properties by 
external magnetic fields. 

We concentrate on anisotropic uniaxial magnetic 
gels^^i^i^. They are manufactured by applying a 
strong external magnetic field during preparation, which 
leads to the formation of oriented straight chain-like ag¬ 
gregates of embedded magnetic particles^”—. After sub¬ 
sequent chemical cross-linking of the embedding polymer 
network, the particle positions get permanently locked^. 
We assume that the magnetic moments carried by the 
particles are free to reorient. First, for diameters up to 
10-15 nm, this applies within the interior of each mag¬ 
netic particle^. Second, this is possible when each parti¬ 
cle as a whole is free to rotate^, e.g., when the polymer 
is not completely cross-linked in the immediate particle 
vicinity^. Another example are yolk-shell particles with 
a magnetic core that can rotate within the shell^i^. If 
reorientations of the magnetic moments are blocked, only 
the first of the two mechanisms described below is active. 

We identify a superelastic stress-strain behavior of uni¬ 
axial magnetic gels by numerically investigating the fol¬ 
lowing model system. Identical spherical colloidal parti¬ 
cles, each carrying a permanent magnetic dipole moment, 
are embedded in a continuous elastic matrix. The elastic 
deformation energy of the matrix is described by a stan¬ 
dard nearly-incompressible Neo-Hookean mode l^^i^^ . We 
tessellate the matrix into sufficiently small tetrahedra by 
Delaunay triangulation^. Each tetrahedron may deform 
afhnely, increasing its elastic energy, from which we ex¬ 
tract restoring forces on its delimiting nodes. Nodes at- 
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tached to surfaces of rigid embedded particles transmit 
forces and torques to these particles. Energy minimiza¬ 
tion with respect to all node and particle positions, as 
well as all particle and dipole orientations, is performed 
(see supplemental material for technical details^). 

Within our model, we study small three-dimensional 
systems, each containing 96 magnetic particles (Fig. [Taj). 
For initialization, we arrange the particles in straight 
linear chain-like aggregates: each chain is one parti¬ 
cle in diameter, but several equi-distanced particles in 
length that are separated by finite gaps filled with elas¬ 
tic material^i^i^-— . The chains are initially aligned 
parallel to each other, but otherwise placed in a random 
non-overlapping way^. Through the presence of the rigid 
inclusions, the elastic modulus increases^ by a factor of 
^ 7. Finally, the magnetic moments are switched on and 
the system is equilibrated, leading to an initial matrix 
deformation (Fig. [Ta|). 

We quasistatically stretch our systems along the chain 
direction. To impose a certain extension, the mesh nodes 
at two opposite system boundaries are displaced into op¬ 
posite directions in small steps. After equilibration dur¬ 
ing each step, the forces on the boundary nodes are mea¬ 
sured. We check the reversibility of the induced total 
deformations by repeated loading and unloading cycles. 
Forces F are measured in units of Fq = ER^ ^ magnetic 
fields B in units of Bq = and magnetic dipole 

moments m in units of mo = R^yJ^E. Here, E is 

the elastic modulus of the matrix, R the particle radius, 
and /io the vacuum permeability. Lq denotes the initial 
total length in stretching direction, AL the (absolute) 
elongation, and AL/Lq the elongational strain. We fix 
the material parameters by setting the Poisson ratio of 
the matrix to = 0.495 and the magnetic moment to 
m = 10 mo- 

Figs. Ilbllcl illustrate the resulting markedly nonlin¬ 
ear stress-strain behavior. First, the force to achieve 
a certain elongation steeply increases with the imposed 
strain. Then a pronounced superelastic nonlinearity fol¬ 
lows. Since our measurements are strain-controlled and 
due to the finite size of our systems, we observe a regime 
of negative slope. A macroscopic sample in this region 
would become inhomogeneous leading to a plateau-like 
superelastic regime^i^ or likewise show negative slope un¬ 
der strain control^. In this area, a slight further in¬ 
crease in applied force induces a huge additional defor¬ 
mation. Remarkably, we can reversibly shift the non¬ 
linearity to smaller strains by a perpendicular external 
magnetic field (Fig.jTbj). High field strengths even switch 
off the non-linearity. Furthermore, we can alter the shape 
of the plateau-like regime by a field in stretching direc¬ 
tion (Fig.[Tc|). At the end of the plateau, the stress-strain 
curve crosses over to a relatively constant intermediate 
slope. 

We found that a combination of two effects allows for 
this adjustable superelastic behavior: a stress-induced 
detachment mechanism of the individual chain-like ag- 




FIG. 1. (Color online) [(^ Snapshot of an equilibrated system 
of chain-like aggregates before stretching. Small bar mag- 


nets in the embedded particles indicate dipole moments, (b) 


Uniaxial stress-strain behavior for stretching in chain direc¬ 
tion, revealing a pronounced superelastic plateau-like nonlin¬ 
earity, and adj ustability by magnetic fields perpendicular to 
the chains. |(c)| Effects of magnetic fields parallel to the chains. 
As in all subsequent figures, solid lines represent loading and 
dotted lines unloading, respectively, highlighting reversibility. 


gregates (Fig. |2|) plus a reorientation mechanism of the 
magnetic moments fFigs. [3l4|) . To illustrate the first one. 
Fig |2a] shows the initialized state of an example chain 
system. Switching on magnetic interactions, the parti¬ 
cles attract each other and “pre-compress” the elastic 
gap material fFig I2b|) . Now the particles are located at 
a small distance from each other, its inverse cube set- 
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FIG. 2. (C olor online) Illustration of the detachment mecha¬ 
nism: |(a)| Cut through an initialized four-particle single-chain 
system, the elastic matrix not yet deformed. |(b)| Equilibrated 
state with pre-compressed gap material. Color maps show 
the local matrix distortion along the longitudinal stretching 
directi on, here illust rated using the so-called “true strain” 
|(c)| First and |(d)| second detachment events of parts 
of the chain where magnetic energy barriers are overcome by 
longitudinally applied stretching forces. |(e)| Equilibrated un¬ 
loaded state in the presence of a strong perpendicular mag¬ 
netic held. All magnetic moments are realigned, thus there 
is no pre-compression and no magnetic energy barrier. m 
Stress-strain curves for various perpendicular held strengths. 
Each spike corresponds to a detachment event where a mag¬ 
netic energy barrier is overcome. An increasing perpendicular 
magnetic held lowers the detachment threshold, until supere¬ 
lasticity is switched off. 


ting the dipolar interaction scale. The initial steep in¬ 
crease of the stress-strain curve (Fig. [JJ) rehects these 
strong magnetic interactions. Once the magnetic barrier 
is overcome, a small further increase in stretching force 
is sufficient to detach part of the chain from the remain¬ 
der (Fig. [2^. Such events suddenly elongate the system 
and lead to spikes in the stress-strain curve (Fig. [JJ. 
They repeatedly occur (Fig. [M|) until all particles have 
detached from each other. In total, a spiky plateau ap¬ 
pears (Fig. [23). It smoothens when averaged over differ¬ 
ent chains (Fig. [1]). This is an intra-chain effect. Inter¬ 
chain interactions are minor for experimentally reported 
particle fraction o^^i^^ . Thus, the finite size of our model 
systems will not affect the significance for the overall ma¬ 
terial behavior. 

Now, the tunability by a perpendicular magnetic field 
becomes clear (Figs. HiEI. Strong perpendicular fields 
align all dipoles in the perpendicular direction. There 
is no pre-compression (Fig. |2e]), hence no magnetic en¬ 
ergy barrier for pulling the particles apart, thus no cor¬ 
responding stress-strain nonlinearity. Smaller magnetic 
fields do not significantly alter the dipole orientations in 
the pre-compressed state due to the strong magnetic in¬ 
teractions at short distances. Yet they affect the thresh¬ 
old for detachment (Fig. jSJ): the moments can rotate 
away from the chain axes when the particle separation 
increases. 

The second effect contributing to the superelastic non¬ 
linearity results from stretching-induced reorientations of 
the magnetic dipole moments. For illustration, we ex¬ 
plain it on a regular cuboid lattice arrangement (Fig.|3|). 
Initially, the edges of the cuboid unit cells are shorter 
along the stretching direction (Fig.jSaj). Thus the dipoles 
align parallel to it. During elongation (Fig. [3b)), these dis¬ 
tances increase. Simultaneously, due to overall volume 
preservation, the system contracts from the sides. This 
decreases the separation perpendicular to the stretching 
direction. When the edge lengths of the distorted unit 
cells become equal in both directions, there is no single 
favored dipole orientation left (Fig. [3b]). Further stretch¬ 
ing makes the dipoles rotate into the plane perpendicular 
to elongation (Fig. [3c]). Thus, during the overall process, 
the magnetic dipoles reorient (“flip”). This flipping ef¬ 
fect is sensitive to the spatial particle arrangement. Yet, 
it likewise appeared in all of our investigated uniaxial 
systems. As becomes obvious from Fig. [Sb] we can iden¬ 
tify it by a vanishing orientational order of the magnetic 
moments (see supplemental material for a quantitative 
evaluation^). A similar orientational analysis for the 
separation vectors between all nearest-neighbor particles 
demonstrates that the mechanism is indeed triggered by 
the distance changes mentioned above^. 

The flipping process is reflected by a steep step in the 
corresponding stress-strain curve (Fig.[ll see supplemen¬ 
tal material for a simplified energetic model including 
larger system sizes^). Also this flipping contribution to 
the superelastic stress-strain behavior can be tuned from 
outside. Through a perpendicular magnetic field, flip- 
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FIG. 3. (C olor o nline) Illustration of the reorientation mechanism: snapshots of a stretched regular cuboid-lattice system (here 
m = 8mo). 1(a) [ Initially, the particles are closest in stretching direction (horizontal), their magnetic moments aligning along this 
direction. Stretching the system increases distances along and (due to volume preservation) decreases distances perpendicular 


Stretching the system further makes a perpendicular direction the preferred one; the orientations of the magnetic moments 
“flip”. 


to the stretching direction. |(b)| When these distances become approximately equal, no particular orientation is preferred, (c) 


ping can be shifted to lower strains or be switched off 
completely, which largely eliminates the corresponding 
nonlinearity in the stress-strain curve (Fig. l4a|) . However, 
also a parallel field has a significant influence: it can post¬ 
pone flipping to larger deformations and smear out the 
connected stress-strain nonlinearity (Fig. I4b|) . In Fig. [T^ 
it is the dip in the stress-strain curve that is mainly con¬ 
nected to the flipping mechanism and can be switched 
off by the parallel field. The steep jumps of Fig. 0] ap¬ 
pear more rounded in Fig. II bllcl because all stress-strain 
curves in Fig. [1] were obtained by averaging over twenty 
characteristic numerical realizations. For larger system 
sizes, the curves would appear still more rounded. 

To amplify the effects and to realize our assumed value 
of m = 10 mo, strong magnetic moments and soft elastic 
matrices are preferred. A possible route is to use parti¬ 
cles made of a material of high remanent magnetization, 
for example NdFeB^ (more than 2 x 10^ A/m). Soft elas¬ 
tic matrices of < 10^ Pa can be made of silicon o^^i^^i^^ 
or polydimethylsiloxane^. The problem is qualitatively 
invariant under rescaling all lengths by a characteristic 
dimension such as the particle radius R. Thus, the par¬ 
ticle size is not a critical factor. Our calculations were 
perfomed for permanent reorient able magnetic moments. 
Yet, the reorientation effect could likewise be observed 
using soft magnetic particles that are magnetized by an 
external magnetic field. Then the reorientation process 
must be performed “manually” by switching the external 
magnetic field direction. 

In conclusion, we have identified a superelastic stress- 
strain behavior of soft uniaxial magnetic gels. These ma¬ 
terials consist of chain-like aggregates of magnetic col¬ 
loidal particles embedded in a soft elastic polymer ma¬ 
trix. Stretching the systems in chain direction reveals 
a significant nonlinearity in the stress-strain curve. In 
this regime, the systems can be strongly deformed with 
hardly any additional load necessary. Two underlying 
stress-induced mechanisms were identified: a detachment 
mechanism of the embedded chain-like aggregates and a 


reorientation mechanism of the magnetic moments. Both 
processes are reversible upon stress release, in analogy to 
the superelastic behavior of shape-memory alloys. As 
an additional benefit, the superelastic properties in the 
present case can be reversibly switched during operation 
by an external magnetic field. 

These nonlinear stress-strain properties open the path¬ 
way to numerous applications. The on-demand tunable 
deformability could be used for easily applicable pack- 
agings or gaskets that become rigid when an external 
magnetic field is switched off. Under pre-stress, external 
magnetic fields can trigger significant deformations, qual¬ 
ifying the materials for the use as soft actuators. Com¬ 
bined with an increased biocompatibility, these concepts 
should be transferable to medical applications such as 
quick wound dressings, artificial muscles, or tunable im¬ 
plants. Finally, combinations of magnetic gels with other 
materials can bestow the tunable superelastic properties 
on the resulting composite. 
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